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We analyze the oscillation signals generated by one extra sterile neutrino. We fully take into 
account the effects of the established oscillations among active neutrinos. We analyze the 
effects in solar, atmospheric, reactor and beam experiments, cosmology and supernova. We 
analyze the impact of the LSND anomaly on cosmology showing, in the full 4 neutrino context, 
that it is not compatible with standard cosmology. We identify the still allowed regions of the 
parameter space, outlining which are the future experiments which can improve the bounds. 

1 Introduction 

The solar and atmospheric neutrino anomalies are today very well explained by oscillations 
among the three active Standard Model (SM) neutrinos. The contribution from a fourth sterile 
neutrino once was an alternative explanation the the observed anomalies, but today it can 
represent only a subleading contribution to the standard scenario of active-only oscillations. 
Many extensions of the Standard Model foresee the existence of fermions which might have a 
mass of the order TeV^/Mpi. A few candidates are, in alphabetic order, axino, branino, dilatino, 
familino, Goldstino, Majorino, modulino, radino. The relevant questions today are then the 
following: how large can be a subdominant sterile neutrino effect and where do we have to 
look for it? Of course the mixing of a sterile neutrino with the active ones affects directly the 
neutrino experiments (solar, atmospheric, reactor and beam). However neutral fermions with 
eV-scale mass, which is the definition of "sterile neutrinos" for us, are typically stable enough 
to give effects in cosmology (Big-Bang Nucleosynthesis, Cosmic Microwave Background, Large 
Scale Structure). Furthermore active-sterile neutrinos mixing would have modified the fluxes 
of active neutrinos coming from the SN1987A supernova with respect to the standard scenario 
(with no sterile neutrinos). 

After presenting our non-standard parametrization of active-sterile mixing (Sec. \^ we will 
summarize the effects of a sterile neutrino on solar and KamLAND experiments (Sec. 13)), on 
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Figure 1: Basic kinds of four neutrino mass spectra. Left: sterile mixing with a flavor eigenstate {v^ in 
the picture). Right: sterile mixing with a mass eigenstate {1^2 in the picture). 



atmospheric, reactor neutrinos, short and long-basehne neutrino beams (Sec. IIJ, on cosmology 
(Sec. El) and on supernova (Sec. The results will be shown in a final set of plots which 
will identify the regions of parameter space which are excluded (in the case of experiments) or 
strongly disfavored (in the case of cosmology and supernovae) . To avoid to show an unreadable 
plot, we show only a few lines and shade the whole region excluded by both solar and atmospheric 
(and reactor and beam) experiments. To better understand the separated bounds we refer to 
the plots o 

2 Parametrization of active-sterile mixing 

Active-sterile neutrinos mixing is usually considered assuming that the initial active neutrino 
|z^a) oscillates with a single large Am^ into a mixed neutrino cos6s\i^'g} + sin0s|fs)- We relax 
these simplifying assumptions and study the more general 4-neutrino context. 

In absence of sterile neutrinos, we denote by U the usual 3x3 mixing matrix that relates 
neutrino flavor eigenstates i^e,fi,T to active neutrino mass eigenstates i^f2 3 f£ = Uni'f {i = 
{1,2,3}, i = {e,/x, r}). In order to parametrize the mixing of an additional sterile neutrino we 
introduce a complex unit 3-versor n which defines the combination of active neutrinos 

n-V = UeVe + n^Vy^ + UrVT = "l^^l + ^21^1 + n^U^ {m = UeiUi) (1) 

to which the sterile neutrino mixes with an angle ^s- With this parameterization the fourth 
mass eigenstate is a superposition of the sterile neutrino neutrino and of the combination 
n • z7: 1/4 = 1/3 cos^s + nivi sin^g. The connection between our parametrization and the usual 
ones can be found in- . In Fig. ^we exemplify the two limiting cases we focus on: mixing with 
a flavor eigenstate (fig. ^) and mixing with a mass eigenstate (fig. E^)- 

From now on we will consider a 3 + 1 scheme, with normal hierarchy and we will set ^13 = 0. 
We have checked that allowing ^13 to assume its maximum allowed value leads to at most minor 
corrections to our analysis. 



3 Sterile neutrinos and solar (and KamLAND) experiments 

The solar neutrino anomaly is today explained by — i^^^r oscillations in the LMA-MSW 
region with mixing parameters tan^ ^sun = 0.41 it 0.05, Am^^^ = (7.1 it 0.6) • 10~^eV^. We 
looked for a subdominant sterile neutrino signal in the experimental data. The work we have 
done is the following. First of all we compute, in the full 4 neutrinos context, the survival 
probability Pee at the detection point for an electron neutrino produced inside the sun. The 
most complicated task is to follow the evolution of the neutrinos inside the sun, which is better 
achieved in terms of the 4x4 neutrino density matrix. Possible non-adiabaticities are considered 
with usual crossing probabilities Pc which however are computed analytically in a non-standard 



way^^. This allows to develop a faster numerical code. Then we fit all the experimental data: the 
SNO CC+NC+ES spectra; the total CC, NC and ES rates measured by SNO with enhanced 
NC sensitivity; the Super-Kamiokande ES spectra; the Gallium rate obtained averaging the 
most recent SAGE, Gallex and GNO data; the Chlorine rate; the KamLAND reactor anti- 
neutrino data. We calculate a as function of the 2 parameters that describe sterile oscillations, 
marginalizing the full with respect to all other sources of uncertainty including the LMA 
parameters Am^^^ and ^sun- This is motivated by the fact that experiments allow only relatively 
minor shifts from the LMA point. 

No statistically significant evidence of a sterile neutrino has been found. The results are 
shown in Fig. |31 where we shaded the region excluded at 99% C.L. Notice that in the same 
plot we shaded also the region excluded at the same C.L. by atmospheric, reactor and beam 
experiments (Sec. For separated plots see^. In^we also outlined a few promising signals. 
The most powerful one would be to discriminate a 2 % shift from the LMA region of the survival 
probability Pee at sub-MeV energies. In fact it is at such energies that LMA oscillations allow 
a vi component in the solar neutrino flux: at higher energies matter effects are so strong that 
the solar fiux is constituted by 1/2 only. The sterile neutrino Vs can mix with vi either directly 
mixing to it or, independently on the mixing, because the two levels cross^. 

4 Sterile neutrinos in atmospheric, reactor and beam experiments 

The atmospheric anomaly is today explained by 1/^ oscillations with parameters sin^ 20atm = 
1.00 ± 0.05, Am^^jn = (2.0 ± 0.4) • lO-^eV^ We looked for subleading effects due to sterile 
neutrinos in the data. We followed the evolution of the 4x4 neutrino density matrix p from the 
production to the detection point. In each medium (air, mantle, core) the evolution is given by 
a diagonal matrix of phases diag exp(— 2i5) where 5i = m'^^.L/AE^. We fitted all most recent 
results and built a global which we marginalized with respect to Am^tm and ^atm- This 
makes the computation much simpler and is motivated by the fact that the experimental data 
allow only minor shifts from the Am^^j^,^atm point. We do not include in our fits the LSND 
data, which have been analyzed separately. In Sec. El we analyze the cosmological impact of the 
LSND anomaly. 

It is useful to distinguish two class of experiments. In the data which are not sensitive to the 
atmospheric anomaly (Chooz, Bugey, CDHS, CCFR, Karmen, Nomad, Chorus) one has 
essentially to handle with vacuum oscillations. Disappearance experiments provide the dominant 
constraints. Instead, the experiments which see the atmospheric anomaly (SK, MACRO and 
K2K) probe sterile neutrinos in a significant but indirect way, essentially by the zenith-angle 
spectra of //-like events with TeV-scale energies. 

It is useful to compare the sensitivity of the these two classes of experiments. Since there are 
no MSW resonances, all these experiments are sensitive only to relatively large sterile mixing, 
^0.1. Sterile mixing with Ve (and with the ui and 1^2 mass eigenstates that contain a sizable 
z/e fraction) is better probed by reactor experiments, although e-like events at SK extend the 
sensitivity down to smaller values of Am?. On the contrary atmospheric experiments give more 
stringent tests of i^s/^'r mixing and of Vs/^^ mixing. 

The results are shown in Fig. IHlwhere we shaded the region excluded at 99% C.L. (as already 
said above, for separated exclusion plots for solar and atmospheric experiments see*^. 

5 Sterile neutrinos and cosmology 

Cosmology can test sterile neutrinos in the following ways. 

Big Bang Nucleosynthesis (BBN). BBN probes the total energy density of the Universe 
at T ^ (0.1 1) MeV. Given a few input parameters (the effective number of thermalized rela- 
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Figure 2: The LSND anomaly interpreted as oscillations of 3+1 neutrinos. Shaded region: suggested 
at 99% C.L. by LSND. Black dotted line: 99% C.L. excluded from other neutrino experiments (mainly Karmen, 
Bugey, SK, CDHS). Continuous red line: A^^ — 3.8 thermalized neutrinos. Dot-dashed orange line: fl^h'^ = 0.01. 

tivistic species, the baryon asymmetry ns/n^ = rj, and possibly the fe/ue lepton asymmetries) 
BBN successfully predicts the abundances of several light nuclei. Today rj is best determined 
within minimal cosmology by CMB data to be^^ r] = (6.15 it 0.25) • lO"^'^. Thus, neglecting 
the lepton asymmetries (which is an excellent approximation unless they are much larger than 
baryon asymmetry) one can use the observations of primordial abundances to test if the number 
of thermalized neutrinos is N^^ = 3 as predicted by the SM. 

In our analysis we fix the active-active oscillation parameters to their experimental values, 
and given the sterile-active mixing parameters we solve the 4x4 neutrino density matrix kinetic 
equations and follow the relevant network of Boltzmann equations in order to compute the 
'^He,D abundances. At the end we convert the computed abundances in an effective number of 
neutrinos. We consider Ni, > 3.8 disfavored by the experimental data^l, while Ni, > 3.2 could 
be tested if the experimental determinations of the ^He,D abundances improve. 

Cosmic Microwave Background (CMB). CMB anisotropics are sensitive to the energy 
content of the Universe at temperatures T ~ eV, since it determines the pattern of fluctuations 
which is measured by WMAP (and other experiments). Neutrinos affect CMB in various ways 

Global fits at the moment imply'^ N^^^ ~ 3 =b 2 somewhat depending on which priors and 
on which data are included in the fit. Future data might start discriminating 3 from 4 neutrinos. 

Large Scale Structure (LSS). Neutrinos can be studied looking at the distribution of 
the galaxies because massive neutrinos move without interacting, making galaxies less clustered. 
Observations can constrain the neutrino relic density Q^h'^ = Tv[m ■ />]/(93.5 eV), where m is 
the 4x4 neutrino mass matrix and p is the 4x4 neutrino density matrix. We approximate the 
present boundl^lwith Qi^h'^ < 0.01. A 0.001 sensitivity in Q^h'^ could be reached in the future. 

For the technical details see^. The results are shown in Fig. |31 The relevant line is the 
red one, which corresponds either to a number of effective neutrinos for BBN equal to 3.8 or 
to a neutrino relic density equal to 0.01, which are both excluded in standard cosmology. The 
bound from BBN dominates at lower 5m^, while the bound from LSS is stronger at higher 6m? 
and small mixing angles. For the i^i/i's mixing the bound is stronger then for the V2fil'^s mixing 
since v\ has the maximal component (see Fig. BBN is more sensible to since it enters 
directly in the weak reactions which interconvert neutrons into protons. For the mixing with 
flavor eigenstates the effect does not vanish for ^rr? — > since the sterile component is present 
in all the mass eigenstates (see Fig. ^). 

We also analyzed the impact of the LSND anomaly on cosmology. Extending previous 
analysis to the full 4 neutrino context, we show our results in Fig. [21 We can conclude that 




the region which explains the LSND anomaly through a sterile neutrino is strongly disfavored 
by Standard Cosmology (BBN primarily) since it would correspond to a number = 4 of 
thermalized neutrinos which seams to be disfavored by ^He data. 

6 Sterile neutrinos in supernovae 

Supernovae arise from the gravitational collapse of stars. The collapse begins when the iron 
core of the star reaches the Chandrasekhar limit. The collapse abruptly stops when nuclear 
densities are reached: the falling material bounces on the surface of the inner core and turns the 
implosion of the core in an explosion of the outer layers. Although there is not yet a definite 
theory of supernova explosion a few key features are quite robust. The beta reactions effectively 
act as a continuous pumping of energy and lepton number from the core matter (which gets 
neutronized) into the neutrinos, that carry them away and will give rise, at the end of the game, 
to a neutron star. Independently on the details which give rise to the supernova explosion, about 
99 % of the gravitational binding energy of the progenitor star (typically 10^^ erg) is released 
through neutrino emission in a few tens of seconds. Despite the low statistics, observations 
from SN1987A show an agreement between these estimates and the observed neutrino flux 
duration and intensity. If one adds channels where the active neutrinos can escape (sterile 
neutrinos in this case, but also for example extra dimensions^ one has to handle with the 
so called "energy constraint": the total gravitational binding energy of the progenitor star is 
pretty well estimated and cannot be drastically overcome. All these considerations allow to use 
supernova to put constraints on the mixing between active and sterile neutrinos. Other bounds 
could be considered. For example in the neutrino-driven picture of supernova explosion the 
active neutrinos flux is essential in order to provide the necessary energy to the outer material 
to successfully explode. However the "energy constraint" seems significantly more robust and it 
is the only one we use. 

We compute the reduction of the flux of electron antineutrinos (which have been measured 
byEl) with respect to the standard case (no sterile neutrinos and active neutrinos mixed with 
solar and atmospheric mixing parameters). The techniques to be used are a generalization of 
the solar ones, with a significant complication due to the fact that in a SN the matter density 
grows from zero up to the nuclear densities, making matter effects very strong. A crucial point 
is the matter density and the electron fraction profiles in the mantle of the star. We used the 
results ofl^. The results are shown in Fig. |31 where we plotted the region which produces a 
suppression of 70% or more of the SN1987A Uf, rate (dashed blue contour). If a new supernova 
like SN1987A is observed in the future much more statistics will be available. Furthermore, the 
theoretical understanding of the supernovae explosion mechanism could significantly improve. 
Then, supernovae are a promising tool to test the physics of sterile neutrinos. 

7 Conclusions 

Neutrino physics is maybe the best field where the interplay between experiments and cosmology 
can be seen. We reviewed the effects of the mixing between active and sterile neutrinos taking 
into account the established active-active mixing. We analyzed effects on solar, atmospheric, 
reactor and beam experiments, cosmology (BBN, CMB and LSS) and supernovae. We found no 
statistically significant evidence of sterile neutrinos and identified the still allowed region of pa- 
rameter space. We also outlined the future promising signals: sub-MeV solar experiments, more 
precise determinations of primordial abundances, new data for CMB, new data and improvement 
of theoretical understanding of supernova. 
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Figure 3; Summary of sterile neutrino effects. The dashed region is excluded at 99% C.L. (2 dof) by solar or 
atmospheric or reactor or short base-line experiments. The red continuous line corresponds either to 3.8 neutrinos 
for BBN (for lower 5m^) or to a neutrino relic density = 0.01 (for higher Sm^), both excluded by standard 
cosmology. The dashed blue line corresponds to a suppression of the SN1987A rate by more than 70%. 
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